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ABSTRACT 

The Wilkinson Microwave Anisotropy Probe mission produces a map of the cosmic microwave 
background radiation over the entire celestial sphere by executing a fast spin and a slow 
precession of its spin axis about the Sun line to obtain a highly interconnected set of 
measurements. The spacecraft attitude is sensed and controlled using an inertial reference unit, 
two star trackers, a digital sun sensor, twelve coarse sun sensors, three reaction wheel 
assemblies, and a propulsion system. Sufficient attitude knowledge is provided to yield 
instrument pointing to a standard deviation (la) of 1.3 arc-minutes per axis. In addition, the 
spacecraft acquires and holds the sunline at initial acquisition and in the event of a failure, and 
slews to the proper orbit adjust orientations and to the proper off-sunline attitude to start the 
compound spin. This paper presents an overview of the design of the attitude control system to 
carry out this mission and presents some early flight experience. 

INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP), the second Medium-Class Explorer 
(MIDEX) mission, was launched on June 30, 2001 as a follow-on to the Cosmic Background 
Explorer (COBE), which made precise measurements of the cosmic microwave background 
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(CMB) that is believed to be a remnant of the Big Bang marking the birth of the universe. 1-4 
WMAP has measured the CMB anisotropy with sensitivity 50 times that of the Differential 
Microwave Radiometer (DMR) instrument on COBE and angular resolution 30 times finer, 
specifically 20 micro Kelvin and 14 arc minutes, respectively, enabling scientists to determine the 
values of key cosmological parameters and to answer questions about the origin of structure in 
the early universe and the fate of the universe. 5, 6 

Since the major error sources in the DMR data arose from COBE’s low Earth orbit, WMAP was 
placed in a Lissajous orbit around the Sun-Earth L 2 Lagrange point to minimize magnetic, 
thermal, and radiation disturbances from the Earth and Sun. WMAP attained its Lissajous orbit 
around L 2 in early October 2001, about 100 days after launch by a Delta II launch vehicle, using 
a lunar gravity assist following three phasing loops, as shown in Figure 1. The trajectory design 
and operations are discussed in detail in Refs. 7-12. 



Figure 1: WMAP Trajectory to L 2 
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The WMAP instrument includes radiometers at five frequencies, passively cooled to about 90°K, 
covering two fields of view (FOVs) 141° apart on the celestial sphere. The WMAP observatory 

executes a fast spin coupled with a slower precession of its spin axis at a constant angle of 22.5° 
from the Sun line to obtain a highly interconnected set of measurements over an annulus between 
87° and 132° from the Sun. The rotation about the spin axis has to be at least an order of 
magnitude faster than the rate of precession of the spin axis; the rates chosen were 2.784 deg/s 
(0.464 rpm) for the spin rate and -0.1 deg/s (1 revolution per hour) for the precession rate. Figure 
2 shows the scan pattern covered by one of the two FOVs in one complete spacecraft precession 
(1 hour), displayed in ecliptic coordinates in which the ecliptic equator runs horizontally across 
the map. The bold circle shows the path for a single spin (2.2 minutes). As the Earth revolves 
around the Sun, this annulus of coverage revolves about the ecliptic pole as shown in Figure 3. 
Thus the entire celestial sphere will be observed once every six months, or eight times in the 
planned mission life of four years. 


North Ealfp tic Pole 



Figure 2: WMAP Scan Pattern 
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Figure 3: WMAP Spin-Scan Concept 

This paper gives an overview of the Attitude Control System (ACS) that acquires and maintains 
the spacecraft orbit, controls the spacecraft angular momentum, provides for safety, in the event 
of an anomaly, and implements the spin-scan observing strategy while minimizing thermal and 
magnetic fluctuations, especially those synchronous with the spin period. More detail can be 
found in Refs. 13-19. 

ACS OVERVIEW 

WMAP uses three right-handed, orthonormal coordinate systems. The Geocentric Inertial frame 
(GCI) is an Earth-centered frame with its xi axis pointing to the vernal equinox, its Zj axis 
pointing to the North Celestial Pole (parallel to the Earth’s spin axis), and yi = Zi x X;. The 
Rotating Sun Referenced frame (RSR) is a spacecraft-centered frame in which the zr axis points 
from the WMAP spacecraft to the Sun, xr is a unit vector in the direction of zr x z 1} and 
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y R = zr x x R . The RSR frame rotates at approximately 1 7day with respect to the GCI frame. The 
body frame is centered at the spacecraft center of mass with Zb axis parallel to the spacecraft 
centerline, directed from the instrument to the solar arrays, ys axis normal to the instrument 
radiator faces, and xb = yp x zb, as shown in Figure 4. As a minimum, all the WMAP ACS 
control modes maintain the spacecraft zb axis within 25° of the Sun, to satisfy thermal and power 
constraints. 



The WMAP attitude is sensed by an Inertial Reference Unit (IRU), two Autonomous Star 
Trackers (ASTs), a Digital Sun Sensor (DSS), and twelve Coarse Sun Sensors (CSSs); it is 
controlled by three Reaction Wheel Assemblies (RWAs) and a propulsion system. Figure 5 
illustrates the WMAP ACS architecture. More detail on the WMAP ACS hardware suite can be 
found in Refs. 16 and 19. 
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Figure 5: Attitude Control System Architecture 

The IRU comprises two Kearfott Two-Axis Rate Assemblies (TARAs), one with input axes 
aligned with the z B and x B axes and the other with input axes aligned with the zb and ye axes. 
This gives redundant rate inputs on the z B axis; the DSS outputs can be differentiated to provide 
rates on the other axes in the event of an IRU failure. 

The boresights of the two Lockheed-Martin ASTs i8 are in the ±y B directions. Each AST tracks 
up to 50 stars simultaneously in its 8.8° square FOV, matches them to stars in an internal star 
catalog, and computes its attitude as a GCI-referenced quaternion with accuracy of 21 arc- 
seconds (1(7) around its boresight axis and 2.3 arc-seconds (lcr) in the other two axes. 
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The Adcole two-axis DSS has two heads, each with 64° square FOV and an accuracy of 1 arc- 
minute (3c) . The centers of the FOVs of the two heads are in the Xb-zb plane at angles of ±29.5° 
from the zs-axis. The CSSs are cosine eyes located in pairs looking outward from the edges of 
the six solar array panels, alternately pointing 36.9° up and 36.9° down from the xb-Yb plane. 

The RWAs are Ithaco Type E wheels each with a momentum storage capacity of 70 Nms. The 
available reaction torque of each wheel is 0.35 Nm, but this is limited to 0.215 Nm by the 
WMAP software to satisfy power constraints. The reaction wheel rotation axes are tilted 60° 
from the -zb axis and uniformly distributed 120° apart in azimuth about this axis. The wheels 
serve the dual function of counterbalancing the body’s spin angular momentum to maintain the 
system momentum (i.e. body plus wheels) near zero while simultaneously applying control 
torques to provide the desired spacecraft attitude. The wheel axis orientations result in all wheel 
speeds being biased away from zero while the spin-scan observing motion is being executed, 
thus avoiding zero-speed crossings that would occur if the wheel spin axes were oriented along 
the spacecraft body frame coordinate axes. 

The propulsion system comprises eight monopropellant hydrazine Reaction Engine Modules 
(REMs) and associated hardware. Each REM generates a maximum thrust of 4.45 N. 

MOMENTUM MANAGEMENT STRATEGY 

The choice of an Li orbit to minimize magnetic, thermal, and radiation disturbances precludes 
the use of magnetic sensing or torquing. Thus, the propulsion system provided for orbit 
maneuvers and stationkeeping is also used to unload accumulated system angular momentum 
after each orbit adjust. These occur several times in the phasing loops but no more than once 
every three months at L 2 to minimize interruptions of science observations. The RWAs can store 
on the order of 70 Nms of angular momentum in non-spinning modes, and a significant fraction 
of this along the zb axis while spinning about this axis. While executing the Observing Mode 


7 



spin-scan, however, the transverse momentum storage capacity (i.e. in the x B -yB plane) is limited 
to 3 Nms, tire amount that can be cycled among the three RWAs at the fast spin rate without 
adversely affecting attitude control. 

Gravity-gradient, atmospheric drag, and outgassing torques are significant in the phasing loops, 
but the accumulated angular momentum of less than 1 Nms per orbit is easily stored until 
removal following orbit maneuvers at apogee or perigee. Solar radiation pressure torque is the 
only sig nifi cant, disturbance torque at L 2 , and the uniform rotation of the spin axis reduces its 
average along the xb and ye axes by more than two orders of magnitude compared to its 
instantaneous value. The only potentially troublesome component is a “pinwheel” torque along 
the zb axis, which might result from an imperfect deployment of the solar array panels. The 
angular momentum is accumulated in inertial space, so it is clear from Figure 3 that the pinwheel 
torque at one point in the orbit leads to a transverse angular momentum one-quarter orbit, or 91 
days, later. This means that any accumulation of angular momentum from the pinwheel torque of 
more than about 0.03 Nms per day would require momentum unloading more frequently than 
desired. 

Pre-flight estimates of the pinwheel torque gave angular momentum accumulation ranging from 
0.0016 to 0.065 Nms per day, depending on the accuracy of deployment of the solar arrays and 
the resulting symmetry of the spacecraft. 14 The worst-case estimate would reach the Observing 
Mode system angular momentum limit of 3 Nms in 46 days, which is highly undesirable. Flight 
data indicates an angular momentum accumulation of about 0.006 Nms per day, which easily 
meets the three-month requirement. In fact, since this is less than 0.03 Nms per day, Figure 3 
shows that the pinwheel torque will begin to unload the accumulated angular momentum on the 
next quarter orbit, so no unloading by the REMs is required at all, in principle. The orbit 
perturbations at L 2 have also been well within requirements, so it has been possible to perform 
stationkeeping and momentum unloading only once every four months, rather than every three 
months. 
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ACS OPERATIONAL MODES 


WMAP has six ACS modes. The Inertial, 
Observing, Delta V, Delta H, and Sun 
Acquisition modes are implemented in the main 
spacecraft (Mongoose V) processor, while the 
Safehold Mode resides in the Attitude Control 
Electronics Remote Services Node (ACE RSN).- 
Figure 6 shows the modes and the transitions 
among them. Anomalous behavior can result in 
autonomous transitions from any other mode to 
Sun Acquisition Mode or Safehold Mode, even 
though these transitions are not shown explicitly. 
Each of the modes of the WMAP ACS will be 
discussed below, including a discussion of the 
sensors, actuators, and control algorithms used in 
that mode. Examples of in-flight performance are 
also provided. 



Figure 6: ACS Mode Transitions 


ACS DESIGN AND PERFORMANCE 
Sun Acquisition Mode 

Sun Acquisition Mode uses the CSS, IRU, and RWAs to acquire and maintain the spacecraft zb 
axis within 25° of the Sun, starting from any initial orientation and with any initial body 
momentum less than [13, 13, 55] Nms. This is a thermally safe and power-positive orientation 
before instrument power-on, and this is the mode entered after separation from the launch 
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vehicle. If the rates at entry to the RWA-based Sun Acquisition Mode exceed those that can be 
handled by this mode, the REM-based Delta H Mode is entered to reduce the rates to an 
acceptable level, after which the spacecraft returns to Sun Acquisition Mode. Transition from 
Sun Acquisition Mode to Inertial Mode can be commanded after the Sun has been acquired. 
Transition to the Mongoose control modes from the ACE Safehold Mode is through Sun 
Acquisition Mode. Normal exit from Sun Acquisition Mode is to either Delta H or Inertial Mode, 
depending on the residual spacecraft spin rate. 

The attitude error signals in Sun Acquisition Mode are calculated by computing the cross product 
of the £un vector computed from CSS measurements with the desired sun vector. The nominal 
desired sun vector points the +zb axis, the solar array normal, directly at the sun. The attitude 
error signals are limited and multiplied by the proportional control gain. The rate error signal is 
the body rate measured by the IRUs. The rate error vector is multiplied by the spacecraft inertia 
and then each axis is multiplied by a derivative control gain. The output of this proportional- 
derivative (PD) control algorithm is three torque commands in the body frame. A body-to-wheel 
reference frame transformation matrix is used to transform these commands to the reaction wheel 
frame. The wheel torque commands are scaled down by a common factor if the largest command 
exceeds the wheel torque capability, so that the torque direction is preserved while the largest 
command is the maximum reaction wheel torque command. 

In pre-launch analysis and testing, Sun Acquisition Mode was found to meet its performance 
requirements for all initial system momentum magnitudes of 55 Nms or less. Because this level 
represented 2a separation rates from the Delta II third stage, a contingency thruster momentum 
unload was unlikely, but would be possibly needed. In addition, there were a few degenerate 
cases (e.g., 180 degrees off of the sun with zero rates) that did not satisfy the requirements, but 
these were deemed too unrealistic to be of concern. 
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Performance on orbit was as expected from the pre-launch testing. After a post-launch 
communications gap, contact was established with WMAP via the Tracking and Data Relay 
Satellite TDRS-W at 21:03, 77 minutes after launch and 10 minutes prior to spacecraft 
separation. At the beginning of this contact, WMAP was still in the Delta II third-stage spin; 
gyro rates were saturated (over 5.3°/sec) and the system momentum magnitude measurement was 
over 90 Nms, as shown in Fig. 7. At 21:13, the yo-yo despin from the third-stage spin occurred, 
and WMAP separated from its booster. Within one telemetry update (16 seconds), the gyros 
desaturated with all axes below l°/second. The measured system momentum dropped in that 
update to approximately 10 Nms. The solar arrays began to deploy 14 seconds after separation. 
The spacecraft reported that the arrays were deployed (all arrays deployed to within 25° of their 
fully deployed state) 8 seconds later. The inset in Figure 7 shows the change in the system 
momentum magnitude measurement as the mass properties of the spacecraft changed during 
array deployment and fuel spin down. The arrays opened fully in four minutes, at which point the 
system momentum magnitude was 7 Nms — well within the maximum level of 55 Nms at which 
Sun Acquisition Mode could acquire the sun. 

Figure 8 shows the CSS and DSS measured sun angles at separation. Note that the sun was out of 
the DSS field of view until 5530 seconds of the plot. The spacecraft acquired the sun within 7 
minutes of separation. The spacecraft was declared separated and safe on the sun at 21:26, 13 
minutes after separation. 
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Figure 7: System Momentum Magnitude at Separation 



Figure 8: Sun Angles at Separation 
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In the design of the Failure Detection and Correction logic for WMAP, Sun Acquisition Mode 
was used in many cases as a first level of safing in the event of a problem. Fortunately, there 
were no on-orbit anomalies that caused an autonomous entrance into Sun Acquisition Mode. 
During the first week of the launch and in-orbit checkout period, Sun Acquisition Mode was 
used as the “base” mode to keep the spacecraft in when other operations were not being 
conducted. Once Observing Mode was checked out, it or Inertial Mode became the base mode. 
The maneuver plan used for each of the calibration maneuvers at apogee and orbit maneuvers at 
perigee commanded the spacecraft into Sun Acquisition immediately after the bum. In all cases, 
the mode satisfied all of its requirements. 

Inertial Mode 

Inertial Mode acts as a staging mode between the other operations of the spacecraft; all 
Mongoose modes other than Delta V have the capability to enter this mode. Inertial Mode can 
either hold the spacecraft in an inertially-fixed orientation or slew the spacecraft between two 
different orientations. This is an RWA- and IRU-based mode, with DSS and AST measurements 
used in an onboard Extended Kalman Filter to update the gyro bias and quaternion error 
estimates. Normal exits from Inertial Mode are by ground command only. The high-level 
requirements for Inertial Mode are to acquire and hold a fixed target quaternion within 
10 minutes and to maintain the spacecraft zb axis within 25° of the sunline at all times. The 
desired orientation is commanded as a desired GCI-to-body quaternion q c 20,2] A slew will be 
executed if this is not close to the current spacecraft orientation. The attitude control is by RWA 
torques, which are computed by a PD controller in terms of attitude and rate errors. 

The attitude errors are expressed as twice the vector part q e of an error quaternion, which is the 
quotient of the commanded quaternion and an estimated quaternion q : 

gc =± 9c®<? _1 (1) 
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In this and the following we use the quaternion product convention of Refs. 21 and 22 rather than 
that of Ref. 20 so that the order of quaternion multiplication is the same as that of the 
corresponding direction cosine matrices. The sign in the equation above is chosen to assign a 
positive sign to the scalar component of q t , which has magnitude close to unity for small 
pointing errors. The small angle approximation is used at all times for the attitude error; this 
introduces inaccuracies for large slews, but does not affect the performance of the controller, and 
as there is not direction preference for the slews, it is not a problem. 

A dynamic attitude limiter is employed in Inertial Mode to enable the spacecraft to meet the 25° 
sun constraint during slews. Before implementing the limiter, the sun line angle reached as high 
as 38° in simulated 45° Inertial Mode slews including high spin errors in the z B axis,. The 
dynamic attitude error limiter calculates an attitude error limit for each axis proportional to the 
error in that axis. This preserves the direction of the resulting slew and prevents the spacecraft 
from violating the sun constraint for any slews with spin angles from 0° to 180°. 

The estimated quaternion is computed by the Kalman filter, with IRU, AST, and DSS 
measurements as input. 22-24 This is similar to Kalman filters employed on several previous 
Goddard Space Flight Center missions, except that the AST produces a measured attitude 
quaternion rather than observed star vectors. This simplifies the onboard software computations 
by removing the burden of star identification and the necessity to carry an onboard star catalog. 18 
The Kalman filter will be discussed in more detail later in this paper. 

The onboard Kalman filter also estimates corrections to the gyro drift rates that are subtracted 
from the IRU-measured body rates to produce the body rate vector. Since the desired rates in 
Inertial Mode are zero, the rate error vector to e is the negative of the body rate vector measured 


The commanded control torque in the body coordinate system is given by: 24-26 

where J is the WMAP moment of inertia tensor, and k A , k p , and k t are the derivative, proportional, 
and integral gains, respectively is nominally zero). The factors of two multiplying the error 
quaternion and its integral reflect the fact that the quaternion errors are half the angle errors in 
the small angle approximation. The controller outputs are the same type of torque commands as 
described in the Sun Acquisition section, and are sent to the RWAs in the same manner. 

In addition to accepting ground-commanded attitudes, Inertial Mode has two "precanned" 
quaternions onboard. One will point the spacecraft +zb axis towards the sun, and the other will 
point the spacecraft +zb axis 22.5° from the sun. The commands have fixed quaternions relative 
to the RSR frame described previously. The RSR frame orientation relative to the GCI frame is 
calculated as a quaternion every minute onboard, and this quaternion is multiplied by the stored 
quaternion at the command time to yield the commanded body attitude relative to the GCI frame. 

Inertial Mode also accepts quaternions from a Command Quaternion Table (CQT), which is a 
table of times and quaternions representing a particular attitude profile for the duration of the 
CQT. The CQT can be loaded with different quaternion sequences, and was used to command 
the gyro calibration slew profile as part of the in-orbit checkout. The CQT was also used to 
command the spacecraft to the proper bum orientations, and to command the spacecraft to hold 
the thrust direction parallel to the velocity vector before, during, and after the maneuvers. 

WMAP was in Inertial Mode for a long period of time during the first week of the mission, 
leading to one surprise: Inertial Mode does not hold the spacecraft in an inertially fixed 
orientation. When the spacecraft was commanded to point at the current sun location (via the sun 
pointing precanned command discussed above), the sunline angle from the spacecraft zb axis was 
expected to change by about one degree per day, consistent with the spacecraft’s motion relative 
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to the sun. Instead, the spacecraft telemetry showed that the spacecraft maintained a sun-pointing 
attitude. The difference between the estimated quaternion and the command quaternion slowly 
climbed at approximately 1 deg/day, and was "reset" to zero whenever we commanded a sun- 
pointing attitude, as expected. However, the telemetered attitude errors and sun angles did not 
reflect this. Inspection of the flight code and ACS simulation code revealed that the portion of 
the simulation used to automatically generate flight software only understands commands and 
attitudes in the RSR frame. When a new inertial (GCI) quaternion is commanded, the flight 
software converts it to an RSR reference frame quaternion before handing it off to the controller. 
If this conversion were done each time the controller was entered (i.e., once per second) or when 
the ephemeris was updated, then the command quaternion within the Inertial Mode controller 
would remain inertially fixed. However, this conversion is only done with the initial command, 
so the controller continues to control to the static RSR frame quaternion. As the ephemeris 
updates and the RSR frame moves, the controller moves with the RSR frame. Thus, commanding 
a sun-pointing attitude results in the Inertial Mode controller continuing to follow the sun. This 
was not discovered in testing primarily because nominal operations never involved rem ainin g in 
Inertial Mode for an entire day. The longest stay in Inertial Mode in a simulation run, build test, 
or acceptance test was an hour or less, so the effect was too small to be seen in the test data. 

Observing Mode 

Observing Mode is an IRU-based mode used for science operations, using the same Kalman 
filter as Inertial Mode. The RWAs provide the angular momentum that preserves a near-zero 
system momentum and the torques that maintain both the 22.5° angle between the spin axis and 
the sunline and the desired sky-scan rates illustrated in Figure 2. The tolerance on both rates is 
5%, and the tolerance on the Sun line angle is 0.25°. Observing Mode differs from Inertial Mode 
in that the commanded quaternion is time-varying, the commanded rates are non- zero, and a 
commanded acceleration and a commanded gyroscopic torque are used in a feedforward loop to 
eliminate attitude hangoff. The attitude error is computed as in Inertial Mode, but the 
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commanded attitude quaternion comes from the desired motion relative to the sunline, rather 
than being fixed. Normal transitions into or out of Observing Mode are via Inertial Mode. 

The commanded inertial quaternion is computed as the product of an RSR-to-body and a GCI-to- 
RSR quaternion: 

#BI = *?BR ® *?RI (4) 

The GCI-to-RSR quaternion is computed onboard from ephemeris models. The desired RSR-to- 
body attitude is more conveniently expressed in terms of 3-1-3 Euler angles 20 ' 21 than in terms of a 
quaternion. The commanded values of the three Euler rates are those that give the desired scan 
rates: 

0 C =-l rev /hour = -0.001745 rad /sec 

S c = 0 (5) 

\j/ c = 0.464 rpm = 0.04859 rad / sec 

The commanded values of the three Euler angles are then given by: 

t 

t 

'o 

d c = 22.5° =0.3927 rad (6) 

t 

Vc = Vo + / VA 

‘o 

where % and \jjr 0 are set by the initial state and 0 O by the desired Sun angle. These 3-1-3 Euler 

angles are converted to the commanded RSR-to-body quaternion q m by the standard 
equations. 20,21 

The rate error vector a> e is the difference between a commanded body rate vector co c and the 
body rate vector a> BI measured by the IRU: 

®e = ®c-®BI (7) 
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The commanded body rate vector tu c is computed from the commanded Euler angles and rates 


by the standard equation: 20,21 
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This commanded rate is modified in the Observing Mode controller so that a large rate command 
change when entering or exiting the mode would not cause the spacecraft to violate the sunline 
constraint. If the rate error is too large, than the rate path saturates the wheel command, and no 
attitude error signal gets into the control torque. The result is an uncontrolled spacecraft attitude 
until the rate errors are small. The initial commanded rate is differenced with the previous 
derived rate, limited, and then integrated to get anew derived rate, as shown in Figure 9. This 
new commanded rate is used to calculate the rate error described above and the feedforward 
terms discussed below. The result is a steady increase in the amount of the initial command rate 
fed into the control system, which generates a spin up (when entering) and a spindown (when 
exiting) Observing Mode. This effect is seen in Figure 10. This shows the initial slew (starting 
just after 1000 seconds on the plot) from the z B axis on the sunline to the 22.5° sunline angle, and 
then the spinup into the compound spin, starting at about 1600 seconds on the plot. 



Figure 9: Commanded Rate Limiting 
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Figure 10: Euler Angles Entering Observing Mode 

One of the feedforward terms is the acceleration the spacecraft needs to follow the commanded 
attitude and rates, and is derived by differentiating the equation for 0 ) c with respect to time. 
Since 0 C , i// c , and 6 c are all constant, the commanded acceleration is: 

cosy/ 0 

= sind c -si 
0 

This acceleration is multiplied by the spacecraft inertia to get a commanded torque. 

The second feedforward torque is the gyroscopic feedforward torque, used to move the system 
momentum around among the reaction wheels as required by its constancy in inertial space: 

T g = w c x H sys (10) 



19 




These feedforward torques are added to the output of the PD controller to generate the complete 
control torque for the RWAs: 

T = /[*,<», + *„( 2q.) + k\ (2q.) +<» e ] +T, (11) 

This control torque is in the body frame, and is treated the same as the control torque created in 
Inertial Mode and Sun Acquisition Mode. 

Initially, the precession rate § in Observing Mode did not meet its 5% accuracy requirement, 
showing a 7% variation at the spin period. This was attributed to an inaccurate value of system 
momentum in the gyroscopic feedforward loop arising from a scale factor error in the RWA 
tachometer signals. Evidence for this was that the magnitude of the system momentum, which 
should be constant, had a 0.4 Nms oscillation at spin period and increased during spin-up by 1.0 
Nms. Comparing a high-fidelity simulation with flight data determined that the oscillation and 
spin-up offset could be removed by a small change in the tachometer scale factors by about 2.5% 
for RWA1 and about 4% for RWA2 and RWA3. After loading these new scale factors, the 
variation of the precession rate was dramatically reduced, as were the spin-period oscillation and 
the spin-up offset in the computed system momentum magnitude shown in Figure 11. A recent 
analysis is in basic agreement with this analysis. 27 

Once all sensors and actuators were calibrated. Observing Mode met the 5% requirement on the 
commanded Euler rates and the 0.25° sun angle control requirement, as shown in the time plots 
in Figure 12 and the nearly perfect circle described by the sun sensor data in the body frame 
displayed in Figure 13. 
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Figure 13: DSS Measurements in Observing Mode 

Delta V Mode 

Delta V Mode, which uses the REMs to adjust the orbit in either the initial phasing loops or for 
L2 stationkeeping, is only entered from Inertial Mode by a command sequence specifying bum 
duration, direction and start time. The desired attitude in terms of either a single quaternion or a 
CQT can be configured either via command or by table load. The desired set of thrusters to be 
used is specified via command. The spacecraft remains in Inertial Mode to slew from the initial 
orientation to the desired attitude for the start of the maneuver, and transitions to Delta V Mode 
at the start time of the requested burn. The only sensors used in Delta V Mode are the IRU and 
RWA tachometers. This mode uses a PD controller to hold the spacecraft to a commanded 
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quaternion attitude while executing the Delta V bum. The output of the controller is transformed 
into thruster firing commands using a pulse width modulator with a minimum pulse width of 
0.04 sec. The desired attitude is held by off-pulsing the primary set of thrusters and on-pulsing 
the others. Normal exit is autonomously to Delta H Mode. 

Anomalous spacecraft motion experienced during some early maneuvers was later attributed to 
torques caused by the evaporation of ice condensed on the back of the solar arrays and sun 
shields. 28 Despite this unexpected phenomenon, maximum pointing errors during the nine Delta 
V maneuvers performed in the first three months of the mission were smaller than predicted (3.7° 
vs. 5.5°), and the imparted velocity increments were accurate to 1%. Less than 15 kg of 
hydrazine propulsion fuel was expended to get to L 2 , about half the amount budgeted for this 
phase of the mission. The 57 kg of fuel remaining for stationkeeping and momentum unloading 
at L 2 will easily support a four-year extended mission. 

Delta H Mode 

Delta H Mode uses the REMs to unload spacecraft system angular momentum, which is 
computed using the RWA tachometers and IRU. It is used primarily upon exit from Delta V 
Mode, but can be commanded from Inertial or Sun Acquisition Mode if necessary, although this 
is not anticipated. The same pulse width modulator is used for Delta H as for Delta V, with the 
exception that all thrusters are operated in an on-pulsing manner for Delta H. If entry was from 
Delta V or Inertial Mode, the ACS autonomously transitions to Inertial Mode after the 
momentum has been reduced to less than 0.3 Nms. If Delta H Mode was entered from Sun 
Acquisition Mode, as discussed above, the autonomous exit upon completion of the momentum 
unloading is back to Sun Acquisition Mode. 
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Safehold Mode 


Safehold Mode is implemented in the ACE, so it can be entered autonomously in the event of a 
Mongoose anomaly. It has two configurations, which differ by the rate information used. The 
first, Safehold/IRU, is a copy of the Sun Acquisition Mode in the Mongoose. The second, 
Safehold/CSS, is a minimum-hardware mode using only the RWAs and CSSs, with rate errors 
being computed by numerically differentiating the position error signals. Because it lacks body z 
rate information from the gyros, Safehold/CSS can tolerate less system momentum than can Sun 
Acquisition or Safehold/IRU Mode. Since the CSSs are insensitive to rotations about the Sun 
line, anti-runaway compensation is applied to prevent the wheels from uncontrolled spinning 
about the satellite’s zs-axis. This is accomplished by applying equal damping torques to the three 
wheels if the sum of their speeds exceeds a pre-set value, thereby suppressing z B -axis rotation 
without applying a net torque in the xb-Yb plane. Exit from either Safehold Mode is by ground 
command only. 

Charged particle flux from extreme solar activity on November 5, 2001 caused a power-on reset 
of the Mongoose processor. The ACS transitioned autonomously to Safehold Mode in the ACE, 
which functioned exactly as designed to keep WMAP safe. The transition to Safehold Mode was 
discovered by operations staff at the next telemetry pass about 12 hours later, and recovery to 
Observing Mode was accomplished within three hours of this discovery. 

ATTITUDE DETERMINATION KALMAN FILTER 

Implementation and Operation 

The WMAP Kalman filter is an Extended Kalman Filter that uses AST and DSS measurements 
to update the IRU-propagated attitude and to update the estimated IRU drift rate. The Kalman 
filter operates in Sun Acquisition or Safehold Modes, but the spacecraft does not use its output 
for attitude knowledge in those modes, which control the spacecraft based on an attitude 
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calculated from the CSS measurements. Sun Acquisition Mode is thus ideal for observing the 
Kalman filter performance before using its output to control the spacecraft. The first time the 
Kalman filter was enabled in Sun Acquisition Mode, parameters such as the measurement noise 
covariance matrix, the initial covariance matrix, and the residual tolerances were increased to 
launch values, since the sensors had not been calibrated in flight yet (Table 1). The larger initial 
parameter values allowed the filter to converge and estimate attitude corrections and IRU bias, 
albeit at less accuracy than possible with a fully calibrated system. 


Table 1: Kalman Filter Noise Parameters 


Noise parameter 

Launch Value 

Nominal Value 

DSS measurement 

0.5 deg 

0.2 deg 

- 

AST measurement 

120 asec 

30 asec 

Gyro rate random walk 

0.01296 asec/s 3/2 

0.001296 asec/s 3/2 

Gyro rate white noise 

7.2 asec/sec 1/2 

0.72 asec/sec 1/2 

Initial covariance matrix diagonal 

(2 0 ) 2 , (40°/hr) 2 

(0.5°) 2 , (207hr) 2 


During each control cycle, the IRU, DSS, and AST data are sampled at slightly different (but 
deterministic) times. The previous control cycle’s estimated attitude quaternion is propagated to 
the current sample cycle time with a first-order propagator 20 using the IRU-measured rates, 
corrected with the previous cycle’s gyro bias estimate. The attitude control is executed based on 
that propagated attitude and measured rate. The DSS and AST data are processed in the sensor 
data processing algorithms, and the resulting measurements (sun vector and attitude quaternion) 
are back-propagated to the IRU sample time with another first-order propagator, using the 
current body rate measurement. The back propagation removes any attitude error that may result 
from the timing differences between the sensor samples, especially at the nominal Observing 
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Mode spin rate. The Kalman filter uses the synchronized measurements to update the propagated 
quaternion, the estimated gyro drift bias, and the propagated covariance matrix. 

Submodes within each control mode are used to determine whether or not the filter will update 
the attitude and gyro bias estimate. In general, if the control mode produces accelerations on the 
body, the Kalman filter is set to propagate the covariance matrix, and the estimated attitude and 
gyro bias are not updated. In Inertial Hold State and Observing Scan/Hold State, the Kalman 
filter is used to update the estimates. 

Attitude Estimate 

The measurement residuals, the difference between the measured quantity and the expected 
(modeled) quantity, are a good indication of how well the Kalman filter is estimating the 
spacecraft attitude. Figures 14 and 15 show a typical one-hour period in Observing Mode at L 2 
for the DSS and the AST, respectively. Much of the DSS residual is dominated by noise from the 
measurement, but all axes show some spin-period frequency oscillations. The short-period 
oscillations are due to inaccuracies in the first-order attitude propagator used onboard (discussed 
below), and the longer period oscillations at the precession period are due to small remaining 
misalignments between the AST and the DSS. 

Figure 15 shows the most striking pattern in the AST xs-axis residuals. Good measurement 
accuracy is indicated by the small amount of noise present (a few arcseconds, as in the zb axis), 
but periodic oscillations with magnitude ±30 arcseconds are quite apparent. There is also a large 
15-arcsecond jump near each peak. These oscillations and jumps are about 90° out of phase with 
the y B -axis pattern. In fact, the large jumps coincide with the zero crossing of the spacecraft body 
rate in each axis. Since the body rates are smooth and continuous, the error is probably in the 
attitude propagated from that rate signal. The first-order numerical integrator used to propagate 
the quaternion may induce errors on the order of 15 arcseconds when the rate signal crosses zero. 
Since the zb component of the angular velocity is always non-zero, there are no large jumps in 
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the z B -ax is residual. So, the spin-period dependent oscillations are most likely more a reflection 
of attitude propagation errors than of actual measurement errors. 



Figure 14: DSS Residuals 
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Figure 15: AST Residuals 
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In addition, the AST boresight is along the spacecraft y B axis, so boresight inaccuracies show up 
as larger, noisier y B -axis AST residuals. The 30 arcsecond magnitude oscillations seen in the x B 
and y B axis residuals of both the DSS and AST residuals is probably due to the inaccuracy of the 
first-order propagators used in processing the data and propagating the attitude. The spikes 
correspond to the zero crossings of the x B and y B body rates. The z B axis has a relatively 
constant, nonzero rate, so propagation errors effects don’t affect the residuals, although all 
residual components show some spin-period frequency oscillation. 

Stray light in the ASTs caused some problems during the phasing loops. Both ASTs lost track 
when the Moon was within a degree or so of the FOV, but only for a few seconds in a spin cycle, 
and only for three spin cycles in any precession cycle. No more than 13 AST readings were lost 
in a precession cycle. There is no Moon, Earth, or Sun interference at L 2 , and the ASTs have 
been routinely tracking 15 to 40 stars in the absence of interference. Based on the definitive 
ground attitude solution, the WMAP attitude determination accuracy (RSS Kalman filter error 

and systematic error) is 1.1 arcseconds (x and z, perpendicular to the AST boresight), and 5.3 
arcseconds ( y , around boresight), all la. This easily meets the pointing knowledge requirement 

of 1.3 arcminutes per axis (la). 

Drift Bias Estimate 

WMAP was launched with default gyro drift bias values of 0 deg/hr in each axis loaded in a 
flight software table. These were used to remove the expected gyro drift bias during IRU sensor 
data processing. The Kalman filter estimates corrections to these table values, and that 
correction, the estimated bias, is removed from the IRU rates as well. After the filter was first 
enabled and converged the day after launch, the estimated drift rates were added to the default 
drift rates and loaded into the IRU table onboard. The filter was reset so that it would only have 
to estimate small corrections to the new table drift rates. The total gyro biases have remained 
stable since launch, indicating good gyro health. In fact, the current estimates for the total gyro 
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biases (default table values plus current Kalman filter estimate), are [-10.9, 1.4, 7.3] deg/hour, 
and are within 0.25 deg/hour of the total estimates that were loaded the day after launch. The 
oscillations seen in the x B axis and y B axis Kalman filter estimates and the larger bias seen in the 
z B axis estimate (shown in Figure 16) suggest a small error in the IRU scale factor, rather than an 
actual variation of the sensor data. 



Figure 16: Kalman Filter Estimates of Gyro Bias 

CONCLUSIONS 

The WMAP spacecraft has to meet stringent requirements on attitude determination and control 
while acquiring the sun, slewing, and performing the compound spin necessary to meet the 
mission science objectives. The attitude control system described in this paper has successfully 
met these demanding requirements far from the Earth where no magnetic field is useful for 
sensing or actuation, and with infrequent telemetry passes. The processor upset on November 5, 
2001 illustrated the importance of having a safemode control capability that is independent of the 
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primary control hardware and software. The flight results show that the WMAP Attitude Control 

System meets or exceeds every requirement. 
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